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ABSTRACT 


The shapes of spectra obtained from measurements of ocean 
turbulence are examined. Kovasznay's equation for high wave 
number turbulence is used to describe the inertial through 
dissipation regions. Power spectra of velocity fluctuation 
recorded in the open sea are used to evaluate constants occur- 
ring in the Kovasznay relation. Semi-empirical forms of high 
wave number energy spectra are then displayed. These computed 
Spectra should define a universal form to which mathematical 


models should adhere. 
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Le 2NERODUCT EON 


A. BACKGROUND 

Turbulence occurring at small scales of the order of me- 
ters to centimeters are features that are of interest to engi- 
neers, oceanographers and meteorologists. Much of this 
interest lies in the transport properties of momentum, energy, 
and heat in turbulence even though the bulk of these fluxes 
take place — a wide range of scales. At larger scales 
Such as encountered in geophysics the turbulence is governed 
by mean flow features and hence is influenced by boundary 
conditions and meteorological and oceanographic variables. 
At small scales, the turbulence becomes more ordered (in a 
stochastic sense) as the influence of the mean flow becomes 
less. Certain properties/features are common from one flow 
to another. At reasonably small scales turbulence iS appar- 
ently governed by "universal" laws. A “universal” law implies 
Pmt turbulence found in the ocean, Paine OY Daboraron, 
in gases, liquids, or plasmas will have the same behavior. 

In 1941, A. N. Kolmogorov set up a mathematical framework 
for small scale turbulence. A statistical mathematician, 


Kolmogorov, nonetheless, based his theory on physical reason- 


mre. AS originally formulated, his theorysapplied to se@aies 
of motion which were isotropic. The theory is known as the 
Kolmogorov theory of isotropic turbulence. This theory pre- 


dicted results which were subsequently experimentally veri- 


med. Specifically Kolmogorov predicted that there was a 





range of scales where energy would fall off with wave number 
to the minus five-thirds power. This range where energy is 
transferred inertially is called the inertial subrange and 
describes scales of turbulence which are separate from the 
mean flow but not small enough to be affected by the viscos- 
mey. in this range 


Gtk) = K' ¢e€ 


243 eee (1) 


Z 
where $(k) is the one-dimensional energy spectrum Ca 


(em), K' is a "universal" constant, € is the rate of dissipa- 


e e e cm 2 i‘ 
tion of turbulent kinetic energy (—_) 
sec sec 





» and k is the 
Wave number (ee): 
cm 

Recent experimental work has shown.that this equation 
holds over a much wider range of scales than Kolmogorov anti- 
Cipated. The law holds at large scales where the turbulence 
Cannot possibly be isotropic. This 1s important to resear- 
chers in geophysical fluids because they are provided with a 
law which can be applied to a wide range of geophysical prob- 
lems. 

Turbulence in the atmosphere and in the ocean almost 
miways exhibits an inertial subrange. JIf K’ is known in the 
inertial subrange then € may be determined by measuring @(k). 

inewevaltiation of ¢.is critical in oceanic turbulence 
and air-sea interaction since if € is known, other parameters 
such as the drag coefficient, C and sensible heat flux may 
be estimated. 


The "universal" law therefore easily allows the computa- 


longer coaetie paramount parameter Of the Ehneery, provided 





that K' ois computed aceuratelyceeln ordem—ito compute kK sac- 
curately, it is necessary to experimentally measure € at 
scales where viscosity is taking energy from turbulence and 
transforming it from kinetic energy to heat energy. It is at 
such scales that most of the turbulent energy is dissipated. 
Viscosity becomes important at very small scales because vis- 
cous forces are effective only over short distances. There- 
fore, in order to measure € semen measurement must be 
made at very Spee scales, scales smaller than encountered 

in the inertial subrange region; that is, beyond the inertial 
Subrange. Such difficult measurements have been made over the 
past few years by Grant et al (1961), Pond et al (1966) and 
Nasmyth (1971), as a result of care in recording, analysis 


and improved sensor technology. 


B. APPROACH 

Theoretically there has not been another success like 
lmaimzorov's inertial subrange prediction. There is no single 
comparable expression for the region beyond the inertial sub- 
range although several forms have been suggested. Because 
of the importance of € to other turbulence investigations, 
further theoretical studies at all scales of turbulence will 
be hampered until this small scale region is understood. 

This thesis, therefore, takes not a theoretical approach 
but an empirical approach. By using recently collected data, 
guideposts will be established for the theoretician by deter- 


mining the constants of a semi-empirical equation of the curve 








desc ri bigme. tic shape of the spectrum at and beyond the iner- 
tial subranee region. Certain parametersewil be wewatnaced 
and these evaluated parameters will form guideposts. Any 
eheory, to be correct, should be able to reproduce (thes alee. 
or the empirical curve. Regression techniques will be used 

to determine the constants of a curve describing the shape of 
the spectrum of velocity fluctuations beyond the inertial sub- 


range region. 


C. KOVASZNAY'S EQUATION 
i» theoretical tpackeround 
There have been several proposals attempting to des- 
cribe the transfer of energy across the spectrum (Hinze, 1959, 
fee 190 ££). In terms of the three-dimensional spectrum E(k), 
one can write 


Z 


(ey 2 ee (2) 


dt 
where T(k) represents, at a given Wave number, the energy 
change due to nonlinear interactions. It is the difference 
between energy increase and energy decrease. When this equa- 
tion is averaged, phase information is lost, and as a result, 
there is not enough information available to solve it. The 
averaging, unfortunately, is absolutely necessary at some 
stage in the analysis. To work with the unaveraged equation 
means that one must know the response of the system to ini- 
tial and subsequent perturbations. Because of Gheunatiiecmon 
turbulent flow, this response and the perturbations are 


essentially unknowable. Therefore, in order to solve the 
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averaged equation, information is put in the form civen cea. 
ly reasonable assumptions about T(k). Such proposals are 


wsually introduced in the integral form (2) namely 


oO 


mk =f 1th) dk (3) 
k! 


mhe physical significance of S(k'), which follows from the 
Uietnition of E(k), is that it is the total flux of enerey 
@Meross the wave number k*. Also following from the definition 
bre T(k) is the result 
Li, m0 (4) 
mmenc inertial swtbrange, because in the inertial subrange 
there is a local steady state. Furthermore, before the dis- 
Sipation region is approached 
S(k”") “=a¥e C5) 
the rate of turbulent energy dissipation. It is in this set- 
mee that we introduce the proposal of Kovasznay. 
2. Kovasznay's Proposal 
The essence of Kovasznay's proposal is that S(k") is 
determined by what goes on near k' only. In other words, he 
proposes that energy transfer is a local phenomenon, and is 
determined entirely by local parameters. From (5) and the 


inertial subrange law, then 


Sy) Sy F 


Sct) 2 EG) aan 2) Seeger @ (6) 


where @ is a constant. Then, from (2) 


e) 


T(k) = 52, stk") |v, = {pany 7 & 


5/2 5 ee 
yk + ZOE k } 


k 
(7) 


1 





The negative sign is needed because k' is a lower limit; at 
the upper limit (k = *), E(k) goes to zero very strongly and 
so gives no contribution to T(k). In the inertial subrange 


ob) 2 (8) 





at 
so that 
T(k) = 2vk* E(k) (9) 
Reom (7) and (9) 
avk7E(k)= - dafaet/2SE 49/24 sp8/? 9/7) (10) 
or simplified, 
Por a solution eee 
E = Age >/3 (12) 
where A is a constant and @ is an unknown function. A me 


is chosen since a k dependence of this nature is desired, at 


least in the inertial subrange. When (12) is substituted in 
oP) 
f= ~ 5/3 iim uae! a 35 


Notice that 


Dare: <= rem: -8/3 
ee ee 5 Agk (14) 
melatkes the first terms of —(13) and (11). For the second 
term of (13) 
oy ais ary En hy 2 
as ak =o) que 
nwa - A a ee JG Cis) 


uF 





or 
1 ig A 
! ae ey eS ne 
Upon integrating 
gil? = me Ri f3 pene rglion eeshiol (17) 


Now at low wave number, from (5) 


Eg = Pile ee (18) 
so that 
a a ee gy 2/2 wf? 
= war! 2 g3/2 (19) 
Thus, at low wave number, 
372 Ee 
9 ee ESTO 
w= 
ox | 
(less en a (20) 
a A 


mech Hs’ the constant of integration occurring in (17). 








Therefore 
. 4/3 
gif? . om f ees lee i) 
a ~ VA 
: 4 3 eles 
Lf we introduce the Kolmogorov wave number he (e/v ) , then 
k Aas 
Te Ey Aen) 3 wis 
ae ee Se ys 
a YA 
or 
ey ee Oe _ k 3 4732 
G = (=) = 1 73 (S ie (22) 
Peplacing (22) in (12) 
E(k) = (87/9 9 fa - aT e Ta i (23) 


a5 








which is Kovasznay's formula. It is to be noted this deriva- 
tion leads to an egquationgwhich dififersc@irom the equation 
@escribed in Hinze (1959). 


The appealing aspect of Kovasznay's result is that it 


fits observations well and is relatively easy to use. How- 
ever, at high wave numbers unrealistic asymptotic regions 
appear. Through the use of well defined spectra, the value 


of a will be evaluated and some precision attached to the 
wave number at which Kovasznay'’s equation is no longer valid. 
ace Spmbreat ion 
Liwawconsistent valtie of G@ is found, then. frome 
it is possible to estimate Ke and the shape of the spectrum 
beyond koe Inertial subrange measurements also allow an 


Sottmate of ke from 


2/3 -5/3 


Ek = Kk" k (24) 
Bemece K' is now reasonably well established. Various reasons 
exist for wanting to know the form of the velocity spectrum 
beyond ko» but a recent one concerns laser propagation in the 


free atmosphere and over waves. The smallest scales of motion 


become very important in this application. 
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i. DESCRIPTIONMOE. DATA 


A. DATA COLLECTION 
Six sets of data were used to obtain the spectral descrip- 
tions of oceanic turbulence. These data were obtained by 
Defense Research Establishment Pacific and analyzed by Nas- 
myth (1970). The data used are from his runs. Sections 1 
through 5 follow closely the description given by Nasmyth 
Org970) in his doctoral thesis. 
1, Area 

These data were the results of two seagoing expedi- 
tions carried out in November 1967 and January-February 1969 
in deep water off the west coast of Vancouver Island, British 
Columbia. Previously the measurements of saa scale turbu- 
lence in the open sea by the Defense Research Establishment 
Pacific, Defense Research Board of Canada were made in 1962 
by a submarine with hull mounted sensors. This limited the 
depth of measurement to about 100 meters. The sensor response 
permitted observations of scales down to Z Or 3 Nal lameters 
defined as "well down into the range of energy dissipation 
Peer acGeus EChiects OF up into the dissipaeron ranee Of wave 
numbers." The data of 1967 and 1969 were obtained by a towed 
body permitting ohservations to depths exceeding 300 meters. 

the m@asturements are initwenced by processes  =cccnmmr ne 
in the upper 330 meters of the ocean and environmental aspects 
of the winter months. The measurements were more or less 


randomly scattered over a ten thousand square kilometer area, 
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part way down the continental slope with ocean depths varying 
from 1000 meters hut usually over 2000 meters. The results 
should not be considered representative of the world oceans, 
but the dynamic features of the area are certainly like those 
in many other Parts Opseme, OcCeCan. 
2. Sensors 

The sensors were mounted on the submerged body which 
towed from the surface ship by a servo controlled winch that 
Gompensated for ship motion to an accuracy of about 30 centi-= 
meters in body depth in normal sea states. The winch could 
frjintain the towed body at any desired depth down to 330 
meters. In the "cycling" mode, the winch may be programmed 
to vary depth at a constant rate over a predetermined depth 
interval (typically 6-30 meters) so that at a constant towing 
speed the body follows a saw tooth path through the water. 

The velocity probe consists of an active element of 
evaporated platinum film on a coutcal glass tip. “Electrivea 
connections to the film are made through leads embedded in 
the glass stem. It functions like a hot wire anemometer. 
Since the maximum dimension of the ag gli ay approximacely G5 
mrltlimeters, a resolution of features as small as 2 milli-— 
Meters £5 possible. A wsetul £requeney response up to veo 
Hz is possible. 

The velocity probe is extremely sensitive to vipra— 
tion and extensive precautions woe taken to eliminate inter- 


herence. “Nonetheless, vibration remains the  Pimitine) + acre 
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in the sensitivity of the velocity system through a midrange 
of the spectrum irom about) iz sboezeei2- 

Another major difficulty in the use of this type of 
probe in an oceanic environment is planktonic fouling. This 
is the reason all successful measurements in offshore waters 
have been made in/’the winter months when the plankton concen- 
tration is the least. When a plankton particle contacts the 
velocity probe, the effect is to momentarily inerease the 
effective thickness of the boundary layer over the film. The 
record then shows a sharp transient indicating an apparent 
Momentary decrease in velocity. 

Timrhempianktoneadhenemonily biuaefly, themewas no 
problem, but if one is stuck, then the probe must be cleaned. 
This is accomplished by a system that spews a high velocity 
reverse flow over the tip of the probe without reducing tow- 
ing speed or otherwise interrupting the experiment. The 
velocity probe outputs are watched and the ynfouling can be 
accomplished in a few seconds by manual control from the 
Ship. 

Under favorable conditions during the winter, a 
plankton particle may be contacted every few seconds, but it 
will only be necessary to wash the probe two or three times 
an hour. Therefore, it is possible to obtain records up to 
30°minutes in length with only minor con@minatien which has 


negligible effect on the spectral content of the velocity 


signal. 
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3 oes 

The possible contamination of the velocity "specraagny 
planktonic fouling has been discussed. [heNcontanin a iono. 
the spectra hy temperature is not so easily determined. The 
Magnitude of possible errors was estimated by indirect 
methods as no direct measurements of the dynamic response of 
the velocity probe to fluctuations in temperature are avail- 
able. Without going into detail, I report that Nasmyth con- 
cludes that the velocity signals are at least predominantly 
real velocity without any major temperature component. One 
qualitative argument is that the power spectrum generated by 
the velocity signals do not exhibit any of the unique charac- 
teristics of the spectra of temperature fluctuations in tur- 
bolent flow [Grant, eter oe Vogel, and Moilliet (1968)] and 
the conclusion is that the signal must be predominantly due 
Bo velocity fluctuations. 

4. Spectra 

The power spectra have been computed from a number of 
samples of velocity fluctuations recorded during the November 
mo7 and January-February 1969 cruises. For satisfactory 
spectral analysis, it is necessary to work with reasonably 
uniform samples of as long a duration as possible. Besides 
the problems already discussed, there are other interferences 
that are frequently encountered in the data collection. Tur- 
pulence OCcurs Im patches Within layers. Since there, ic mu0 
way of predicting in advance what the signal levels in the 


next patch will be, many good samples are lost because of 
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inappropriate gain settings. When probe washes are necessary 
it may ruin an otherwise useful Section of LecoOrd. =m ea. 
preciable change in ship's speed during recording will dis- 
ort the resulting spectrum. The point 1S hate it somes 4 
able to get suitable samples of any length, the sample is 
valuable. 

During “Novenber 1967, 32 hours of recording durin. 
two weeks at sea were accomplished. noe January-February 
Jog, 35 ates of recording during three weeks at sea were 
made. From these raw data samples, there are 20 or 25 samples 
varying in length from 10 seconds to two minutes from which 
meaningful spectra were derived. In the cycling mode sample, 
times are short because of changing turbulence layers. In 
the constant depth mode sample, lengths are longer and 
several spectra of 52 seconds duration or longer were obtained. 
Wherever possible, noise is subtracted from the original 
Beeetrum resulting in a noise-free spectrum. 

a. Re COmarm 2 

Components of circuitry which must be located close 
to the sensors are enclosed in pressure cases in the towed 
body. From there, information from all sensors is fed up the 
armored multi-conductor towing cable for further processing 
and recording on the ship. The lower 100 meters of cable are 
faired to reduced vibration. 

On board ship the signal from each thermistor is 
divided into two channels. OQne is amplified and recorded 


direcuby amd =the other is electronically differentia teamwasem 


go 





respect to time befor@® recording..W&his methodic >a sac. 
ceptable sipnal-noise ratio at the higher -=frequencteceoehcr— 
wise unobtainable on a single channel. The signal from the 
velocity probe is divided into three channels™ior "sna 
reasons. The first channel or "A" channel is recorded broad- 
band with no filtering. "B" channel is effectively differ- 
Patiated up to about OMHz, an@® cits of f St @abiewt 20 Hz. ce 
channel differentiates up to 1000 Hz and then cuts off. The 
Signal from the temperature probe is treated similarly except 
that there is a low-frequency cutoff on "A" at about 0.01 Hz 
amd) "B'' extends up to about 100 Hz. 

Signal channels from all sensors together with coded 
timing signals and a voice commentary are recorded on two 
seven channel magnetic tapes with multiplexing as required. 
At the same time, selected channels are monitored on oscillo- 
scopes as an aid in selecting appropriate gain settings and 
12 channels are recorded on two paper charts giving a perman- 
ent visual record. This served as a planning aid during the 


experiment and for later use during analysis. 


a DATA USED 
Data sets 4, 5, and 6 consisted of 45 to 47 data points; 
maereas, data sets 7, 8 and 9 consisted of 174 to 176 data 


points. These six sets of data were plotted on log-log graph 
POO vs log ce )> (Fieures 
(ev ) Ss 

1-6). The data are presented in this form to make comparison 


paper in normalized form as log 


of spectra easier. Normalized plots should lie on top of each 


Z0 





other. The plots covered 8 to 11 decades along the $(k), or 
energy axis and 4 to 6 decades along the k, or wave number 
axis. This remarkably large range of values makes them 
ideally suited) for examining "details of the high vwaveumumber 
region of the turbulent energy spectrum and testing Kovasznay's 
equation. 

The data points clearly indicated the expected -5/3 slope 
mmertial subrange region and the “tailine off" in the Pow 
energy, high wave number region beyond the inertial subrange. 


Several points were discarded by this author because of ob- 


Vious inconsistencies, Only runs 4, 5, and 6 had no points 
discarded. Run 7 had eleven points discarded; 8 had ten 
points discarded and run 9 had nineteen points. Because of 


the similarity of these data sets, it was considered necessary 


Baily to use data set 4 for analysis. 
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TLL. ANALYSIS AND APPROACHES 


me APPROACHES 
1. Litiear Reemess ion 
Since a direct supecei1tution of @ Lor (23a ued — 
cult, the equation was expressed in a form which was amenable 
to regression techniques. Despite several attempts, consis- 
tent results could not be obtained and this method of approach 
to the problem was discarded. 
2, “lierenak Subaange 
An alternate approach was to consider only the ISR. 
In the inertial subrange the one-dimensional energy spectrum 
and the three-dimensional energy spectrum have the same power 
law behavior. The Kovasznay equation then mate be reduced to 


6 (k) : 7 2/3 -2/3 2/3 (25) 


[where @(k) is the one-dimensional energy equation] since 
clearly these terms will dominate. This reduced form allows 
observed @(k) and k to be used in this region to evaluate a. 


Since 


AS 2/3 


23 


O(k) = K' € ac, 


Poem it is natural to associate a with K'. Since K' * 
0.5, then an q@ value of 1.59 can be expected. 
J. Exponential Curve 
Beyond the inertial subrange, @(k) falls off expon- 


entially with k. A desirable equation would be one that 


collapsed to Equation (26) in the inertial subrange but 


Jip 


generated an exponential decay at high wave numbers. Further, 
such an equation should be expressed in terms of meaningful 
physical quantities. This is the goal of the theoretician if 
mot this thesis. The object here is simply to provide 
empirical clues. Consequently, an exponential regression 


analysis was attempted. 


B. ANALYSIS 
1. Linear Regression 
Equation (23) was taken and rewritten in the form 


ee sg ee SPN a Sie 


S 


$8 (k) =a +.25 
(27) 
eamplified into an equation suitable for regression techniques 


this becomes 


ie + A go t/3 + alk (28) 


Y= Ay 2 3 


where Ay» A, and A. are given below. 


The regression analysis program appeared to be very 
momcttiarve to badedata. Accerdane tort hesKovascznay equation, 
A, and A. should be > 0 and A, < 0. Furthermore, [Aa 


> eS Mer such *censistency Magiehoundm@boer thensixecetemot 


| > [4,| 
data when using the Regre Program. Of the six sets, there 
Bec Jlleaet three pains. Each paln corresponded to spectra 
from a single time series of turbulence data. The difference 
SSerelat Ones set wok wre pa i. contained four times as many 
Points as the other. “Even in tite watved dabatcets . calms 


proved inconsistent. Further, the poor regeneration of the 


23 





curve using the Ay> Ay» and A. computed by Repweseausedae aie 
method to be discarded. 
2. Subs Eve ue rom 


Equation (28) was used where 


aes ae i. 250-3 
Bp SUD se "Ie Be 7 eS 
Cak_) ak 
S Ss 
A direct substitution of a = 1.6 was used. The vatwes 
of Ay» Ay» and A, were computed. The values are presented 


in Table I. These values give a regeneration of the curve 
and is an alternative but simpler way of expressing Kovasz- 
nay's equation. A comparison of the curve generated with 
these values of Ay> Ags and A. and actual data is shown in 


3 


Figure /. The agreement is good in the inertial subrange. 


3. Gompertz Curve 

The Gompertz curve y = 2 (atBp”) was used because it 
was readily available. It is unsatisfactory, however, be- 
cause even if solved and a good curve fit obtained, the quan- 
tities a, 8, p, and x do not lend themselves to physical 
interpretation. There are no dimensions on the right hand 
Side. Another curve should be used but the form is not imme- 


diately obvious and is perhaps worthy of a separate research 


investigation. 
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IY. CONCLUSIONS 


Since approaches III-A-1 and ITI-A-3 were unsuccessful, 
poly the result of Approach LIl=A-2 will be discussed 

As mentioned in III-A-2, an a value of 1.59 was expected. 
A spread of a from 1.8 to more than 6.0 was obtained from the 
real data in the manner described below. 

A plot was made of a@ versus wave number Gueene Sue! Wek 
minimum occurs in the vicinity of k = 0.5. Beyond this wave 
number a, values continue to increase with increasing wave 
number. 

Various values of a obtained using the second method 
above were substituted in the full Kovasznay equation and the 
resulting curves compared with actual data (Figure 9). Of 


the several’ curves generated, low values of a generated better 


fits, the best fit occurring for the lowest a, namely a = 1.8. 
This is also closest to the expected value of a = 1.6. The 
fit of curves with a = 1.6 and o = 1.8 were very similar 


although the lower a value curve dropped off more steeply. 
Another feature of interest is that in the inertial subrange, 
O(k) values are lower for higher a (since a is raised to a 
negative power) but beyond the inertial subrange (where other 
terms dominate) the opposite is true. The crossover point 
occurs a¥round k= 20.) Une "sigpmiticanecmoe «Elvis tb ail men 


not clea. 


Zo 





The limit of the Kovasznay result appears to be at a wave 
mumbexy k ~ 50. Beyond thijsewaye number, the yaluessotuencus, 
increase with increasing wave number, a phenomenon that was 
expected, This phememenon obyiousty is inconsis teneweeeee 
served results and indeed would be catastrophic theoretically. 
Lowentere, practically zt 15 not important since present tec 
nology does not allow measurements to be made at such high 
Wave numbers. 

in summary, as Q@ increases, the fitted curve approximates 
Ehe actual curve with decreasing accuracy. The best @ occurs 
at the minimum in the a versus k (Figure 8) and is approxi- 
mately where k < 0.5. The reconstructed curve agrees well up 
meme =~ 2.0. For k > 2.0, @(k) in actual data falls .off more 
rapidly than $(k) computed. Good agreement in the ISR is 
expected and observed. The dissipation range was not expected 
to agree. Figure 9 displays the quantitative results showing 


existing agreement. 
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Figure 1 
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Figure 5 
Velocity Spectrum 
For Data Set 8 





42 





e Figure 6 
Velocity Spectrum 
For Data Set 9 
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